The KATRIN Pre- Spectrometer at reduced Filter 
Energy 



°: 

M. Prall 1 , P. Renschler 2 , F. Gliick 34 , A. Beglarian 3 , 

H. Bichsel 5 , L. Bornschein 3 , Z. Chaoui 6 , G. Drexlin 3 , 

F. Frankle 3 , S. Gorhardt 3 , S. Mertens 3 , M. Steidl 2 , 

CN ■ Th. Thiimmler 3 , S. Wiistling 3 , C. Weinheimer 1 , 

S. Zadorozhny 7 

^j ■ ^Universitat Minister, Institut fur Kernphysik, Wilhelm-Klcmm-Strafie 9, 

^ ! 48149 Miinster, Germany 

I ' 2 Universitat Karlsruhe (TH), Institut fur Experimentelle Kernphysik, 

^ ■ Postfach 6980, 76128 Karlsruhe, Germany 

3 Karlsruher Institut fiir Technologie, Institut fur Kernphysik, Postfach 3640, 
Y\ 76021 Karlsruhe, Germany 

•rH ! 4 KFKI, RMKI, H-1525 Budapest, POB 49, Hungary 

5 CENPA University of Washington P.O. Box 354290 Seattle, WA 98195-4290 
J^, ■ 6 Laboratory of Optoelectronics and Devices, Faculty of Science, 

i *^-f , University of Setif, Algeria 

7 Institute for Nuclear Research of Russian Academy of Sciences, Moscow, Russia 

^. E-mail: matthias.prall@uni-muenster.de 

^": 

^-j- ■ Abstract. The KArlsruhe TRItium Neutrino experiment, KATRIN, will 

C\) . determine the mass of the electron neutrino with a sensitivity of 0.2 eV (90% C.L.) 

^C) \ via a measurement of the /3-spectrum of gaseous tritium near its endpoint of Eq — 

18.57 keV. An ultra-low background of about b = 10 mHz is among the requirements 

to reach this sensitivity. In the KATRIN main beam-line two spectrometers of MAC-E 

. . ' filter type are used in a tandem configuration. This setup, however, produces a Penning 

. J^ ■ trap which could lead to increased background. We have performed test measurements 

S^ | showing that the filter energy of the pre-spectrometer can be reduced by several keV 

H ■ in order to diminish this trap. These measurements were analyzed with the help of 

a complex computer simulation, modeling multiple electron reflections both from the 
detector and the photoelectric electron source used in our test setup. 
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1. The KATRIN Experiment 

The KATRIN experiment pQ will determine the mass of the electron antineutrindjj 

mi, e = JJ2 mj | C/gj | 2 via a high precision measurement of the /3-decay kinematics 
at the endpoint Eq = 18.57 keV of the /3-spectrum of tritium with a sensitivity 
of 0.2 eV (90% C.L.). The latest upper limits obtained with this model- 
independent method and the isotope tritium as /3-emitter are from the experiments 
at Mainz: m^ < 2.3 eV (95% C.L.) [2j and Troitsk: m£ ff > < 2.05 eV (95% C.L.) [3]. 
Figure [1] shows a schematic overview over the 70 m long KATRIN setup: T 2 -gas with 
an activity of 10 11 Bq is recirculated in the so-called windowless gaseous tritium source 
(WGTS). /3-decay electrons are guided by a magnetic field towards the spectrometers. 
Both the pre- and main spectrometer are of the MAC-E filter type [3]. 
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Figure 1. Overview over the KATRIN experiment with the potentials and magnetic 
field strengths given in the KATRIN design report [T]. With these values, the setup 
contains a Penning trap for electrons between the pre- and main-spectrometer at 
I?sol = 4.5 T. A change of the pre-spectrometer potential from C/pg = —18.27 kV by 
several kV towards smaller absolute values will diminish this trap. 



Electrons are guided through these spectrometers along the magnetic field lines 
and are decelerated by an electric filter potential Ut in each of the spectrometers. Ut 
is the high voltage with negative polarity applied to the spectrometer, or more exactly, 
to its inner electrode system pp. In KATRIN, we use the independent filter potentials 
£/f = £/pg for the pre- and Uf = Uy^^ for the main-spectrometer. Only those elec- 
trons with an energy E larger than the filter energy qU^ (see Section |2]) are transmitted 
through a spectrometer and are reaccelerated to their original energy. Here q — — e is 

J KATRIN will not be able to resolve the different neutrino mass eigenstates m- v but will determine a 

weighted average m\ e / of the neutrino mass states m- x according to their mixing C/ e j with the electron 
neutrino. 
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the negative electron charge, and the filter energy qU^ is the maximum potential energy 
of an electron in the spectrometers. The pre-spectrometer (PS) with an energy resolu- 
tion of Ai?pg ps 100 eV is the first filter for the /3-decay electrons. It reduces the flux 
of /3-decay electrons into the main-spectrometer (MS), lowering the rate of background 
electrons created in collisions with residual gas molecules pQ. Having an energy reso- 
lution of AEy^ = 0.93 eV, the MS scans the last 30 eV of the T 2 /3-spectrum which 
contain the information on the neutrino mass. Finally, the electrons transmitted by the 
MS are counted by a 148 pixel PIN diode with an energy resolution of AE^^ ps 1 keV. 

The motivation for our investigations is the following: Inside the MS, the /3-decay elec- 
trons can start multi-step processes leading to free electrons. These can be accelerated 
towards the detector and to energies around 18 keV by the MS potential. The energy 
resolution of KATRIN's detector is about AEu e + ps 1 keV. Therefore, these electrons 
cannot be distinguished from signal electrons produced by the tritium /3-decay. Thus, 
the background can rise above KATRIN requirement of b = 10 mHz [lj. Therefore, 
the flux of /3-decay electrons into the MS should be kept low. The flux of /3-electrons 
can be minimized by keeping the filter potentials of MS and PS relatively close (e.g. 
U FS = -18.27 kV and J7 M g « -18.57 kV). However, using the B-fields (5 gol = 4.5 T, 
Fig. [TJ mentioned in the KATRIN design report jl] the region between the two spec- 
trometers is a Penning trap for electrons. By multi step processes this trap can lead to 
increased background as well [5]. A reduction of the pre-spectrometer potential from 
C/pc = —18.27 kV by several kV towards smaller absolute values will diminish this trap. 
The optimum value of t/pg which minimizes the background has to be determined ex- 
perimentally. If qU-pg is reduced by several keV however, the /3-electrons with energies 
close to Eq = 18.57 keV will retain a surplus energy E snT = Eq — qUpo in the order 
of several keV inside the PS. Thus, electrons have higher speed and may no longer be 
guided by the magnetic field. This behaviour, leading to transmission losses, was al- 
ready observed in the MAC-E filter of the Mainz Neutrino Mass Experiment [61 UJ. For 
KATRIN, 100% transmission above the PS filter energy qUpg is required. We present 
two main results in this publication: 

(i) The requirement of 100% transmission at reduced PS filter energy ^t/pg is fulfilled. 
If it should turn out that the Penning trap between the two spectrometers cannot 
be suppressed by other means [5] , the PS filter energy can be reduced by many keV 
in order to overcome this problem. 

(ii) The KATRIN collaboration is able to model the electron transport and the electron 
backscattering at the detector with high precision in agreement with experimental 
data. This allows detailed investigations of the experimental setup. 

This publication is organized as follows: In section |2] we review the operation principle 
of the MAC-E filter, Section [3] presents our experimental setup; our measurements and 
their analysis via custom simulation tools are presented in Sections H] to El Section [7J 
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discusses the background and adiabaticity at reduced PS filter energy in detail. Finally, 
our findings are summarized in Section [HJ 

2. The MAC-E Filter Technique 

In this Section, we explain the principle of a MAC-E filter [4] under standard conditions, 
i.e. when the filter energy qU^ is close to (a few 10 eV) the energy E of the incoming 
electrons. 
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Figure 2. Left: Working principle of a MAC-E filter. Right: Ideal transmission 
function of the pre-spectrometer for an isotropic source of monoenergetic electrons with 
Eq = 18.57 keV, -B ana pg = 0.016 T and B BQ y — 4.5 T. The transmission probability 
of the pre-spectrometer (of a MAC-E filter) starts to become non-zero as soon as the 
energy E of the incoming electron is larger than the filter energy qUpg. Both the 
tank and electrodes of the KATRIN pre- and main-spectrometer are on negative high 
voltage. Detailed explanations can be found in the main text. 



The principle of a MAC-E filter is illustrated in Fig. |2) Two identical solenoids 
provide a guiding magnetic field. Ignoring drift motions which appear as higher 
order corrections [5], the /9-electrons enter the MAC-E filter and follow the guiding 
magnetic field lines along helix-like trajectories resulting from the cyclotron motion. 
This statement is true if the relative changes of the electric and magnetic field strength 
within a cyclotron length Zcyc are small [8]: 
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Here, v\\ is the electron velocity parallel to the guiding magnetic field line, u cyc the 
cyclotron frequency, tuq the electron mass, q = — e the negative electron charge and 7 
the relativistic factor. 
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If Eq. [TJ holds, there is a conserved adiabatic invariant 7/1, where \i denotes the 
orbital magnetic moment of the electron (see Section 12.5 in [9]) 

l±l.«L = ^EL_ (2) 

p 2 7 B 2m el B y ' 

E± = E^ n -sm 2 (9) is the fraction of the kinetic energy E^ n which can be attributed 
to the motion around the guiding B-field line. 9 is the polar angle between the guiding 
B-field line and the electron momentum vector. In the following, we will use the symbol 
ip for the corresponding azimuthal angle. E\\ = Ekw. ■ cos 2 (9) is the fraction of the 
kinetic energy connected to the forward motion of the electron. Only Eh is analyzed by 
the MAC-E filter. p± denotes the fraction of the electron momentum perpendicular to 
the guiding B-field line. 

For KATRIN, the maximum electron energy is Eq = 18.57 keV, thus one has 
7 < 1.04. Therefore, \x is a good approximation for the conserved quantity, especially 
when electrons are slowed down by the electric field in the MAC-E filter: 

E± E\,:„ • sin 2 9 , , 

fj, ps -± = Km w const. (3) 

B B 

From Eq. [3] it is clear that the polar angle 9 and E± are completely determined by 
the B-field and the kinetic energy of the electron. As B decreases towards the analyzing 
plane of the MAC-E filter, E± is minimized, providing the good energy resolution of 
the MAC-E filter. Electrons are guaranteed to be transmitted along the guiding B- 
field line if their initial energy is large enough to overcome the spectrometer potential. 
The transmission probability T(E, qUf) of a MAC-E filter is derived from Eq. [3] by 
integrating over all electrons which fulfill < E\\ = (E — qUt) — E±_ in the central plane 
of the MAC-E filter. These electrons start at ground potential. Isotropically emitted 
electrons with starting energy E and 9 from 0° to 90° have to be considered. Eq. [2] is 
used to transform between E± at the origin of the electrons with B s and the analyzing 
plane, with the minimum B-field strength £?ana : 



E-qUf B s B s 

T(E,qUf) = l-Jl ^--ET- ^ qU f <E<qU r - |— . (4) 

V ft -t>ana -ds — -Dana 

Thus, the transmission probability T(E, qU£) only depends on the magnetic field 
strengths, on the energy E of the incoming electron and on qUf, the filter energy. Below 
the interval specified in Eq. HJ the transmission probability is zero, above this interval 
it is unity (Fig. EJ). The transmission function describes an energy high-pass filter, only 
electrons with an energy E above the filter potential qll^ are transmitted. The energy 
resolution of the MAC-E filter AEt is equal to the maximum E± of an electron in its 
analyzing plane. E± in the analyzing plane is maximal if the polar angle 9 is equal to 
90° in the entry-side magnet. The conservation of \x allows to compute the resolution 
AEf. 

, E± -Bana /c -n 

const, ps fi = — =^- AEc = E ■ —- — . (5) 

B Ba 
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Considering electrons starting in the entry-side solenoid of the PS, one has to insert 
B s = £? go i = 4.5 T and -B ana pg = 0.016 T. The formula gives an energy resolution 
A£pg = 64 eV if the PS is used with E = 18 keV electrons (cf. right side of Fig. EJ. 
For a MAC-E filter and in the adiabatic approximation [8], the Lorentz force 

F = q(E + vx B). (6) 

does not only result in a helix-like motion around the guiding magnetic field line but 
also in an azimuthal magnetron drift around the spectrometer symmetry axis [4]. 
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Figure 3. Pictorial explanation 
of the Grad BxB-dift. Left: 
The orbit of a charged particle in a 
uniform B-field. Right: The orbit of 
a charged particle in a non-uniform 
B-field. 



Electrons which pass the MAC-E filter on its symmetry axis gyrate around the 
central magnetic field line in a helical cyclotron motion. Since the magnetic field is 
axially symmetric and the cyclotron radius is changing slowly, the electrons are in a 
quasi-constant magnetic field. For off-axis electrons, however, the magnetic field is 
asymmetric during a cyclotron motion. The radius of curvature of the electron trajectory 
is smaller in a stronger B-field. Therefore, the gradient \7B results in an azimuthal 
B x V-B drift of the guiding center along a circle with constant magnetic field (Fig. |3]): 
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(7) 



The lower qU-pg is in our investigations, the larger are E± and E\\ and the azimuthal 
drift (j7|). The following section will show that our PS data can only be understood if 
the azimuthal drift ((7|) is taken into account. 



3. The Pre-Spectrometer Test Setup 

This section describes the experimental setup used to measure the PS transmission with 
E — 18 keV electrons and a PS filter energy qUp^ down to 1 keV, so that electrons 
retain a surplus energy -E'sur = E — qUp^ of up to 17 keV in the PS. 

The B-field of the PS test setup (Fig. Hj) is generated by two solenoids (B). The 
inner electrode system of the PS tank (G) consists of four parts: The ground electrodes 
(C) define the potential at the entry and exit of the PS, the shielding electrodes (D) 
were introduced to avoid a Penning trap leading to background [10] , electrodes (E) are 
conical metal shields, the central part (F) is a wire electrode. The tank (G) and the 
shielding electrodes (D) are electrically connected. Both the tank and the electrodes 
are on negative high voltage. There is a longitudinal gap splitting electrodes (E) and 
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Figure 4. The pre-spectrometer test setup: A) Photoelectric electron source (E- 
GUN), B) solenoid, C) ground electrode, D) shielding electrode, E) full metal cone 
electrode, F) wire electrode, G) spectrometer vessel, H) 64-pixel silicon PIN diode. 
There is a longitudinal gap splitting electrodes E) and F) into a left and right half. 
For our measurements, the left and right dipole half were electrically connected and 
therefore supplied with identical voltages. The dark blue arrow illustrates an electron 
trajectory. 



(F) into a left and right half. For our measurements, the left and right dipole half 
were electrically connected and therefore supplied with identical voltages. Previous 
measurements ^0\ with this setup and a pressure of 10~ 10 mbar inside the PS resulted 
in an average background rate of 17 ± 0.4 mHz in the energy window from 15 to 21 keV 
over the whole detector. For these, the tank (G) and shielding electrode (D) were kept at 
[/tank = —18 kV, the inner electrodes (E) and (F) were put on ^ e j ec trode = — 18.5 kV. 
As this configuration does not produce any background related to particles stored in 
Penning traps, the same potential difference ^electrode — ^tank = — ^.5 kV between 
the tank (G) and the electrode system (E,F) was used in our measurements. The 
potential inside the PS tank is a mixture of the electrode and tank potential, one has 
^PS = a ' ^tank "*" ^ ' ^electrode' ^ e ^' ^ ne cons t an ts are a = and b = 1 in good 
approximation. £7pg and ^electrode never differ by more than a few tens of V. This 
effect depends on the electron trajectory in the PS and is accounted for in our simulations 
(Section E]). For these, we have computed the actual electric field inside the PS using the 
methods described in [lUEG]. I n the following we do not distinguish between gf/pg and 
^electrode * n ^ e ^ ex ^ as their difference is negligible at keV surplus energies. In each 
of our measurement series, the tank voltage £4 an k was var i e d from -0.5 kV to about 
-17.5 kV, so that the PS filter energy qUp^ ~ <7^ e i ec trode var i e d from about 1 keV to 
about 18 keV. 

A photoelectric electron source (E-GUN) (Fig. [5]) mounted at the entry of the PS 
test setup (Fig. H|) was used to generate electrons with an energy of E — 18 keV for the 
measurements: A deuterium lamp (f ) generates UV light with wavelengths in the range 
185 nm < A < 400 nm (6.7 eV > hv > 3.1 eV). The light shines through a sapphire 
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window (e) and a hollow ceramic insulator (d). The UV photons finally produce free 
electrons via the photoelectric effect in a thin gold layer (« 35 /zg/cm 2 ) on a gold plated 
quartz tip (b) sitting in a metal housing (c). The quartz tip is transparent to light 
of the wavelength 150 nm < A < 4000 nm, thus not cutting into the UV-spectrum 
of the deuterium lamp. The work function of gold is 4.83 ± 0.02 eV [13]. Therefore, 
only electrons with excess energies of up to -Emax = 6.7 eV — 4.83 eV = 1.87 eV can 
be released. The gold plated tip is supplied with a voltage of U^ = — 18 kV. The 
photoelectrons are finally accelerated to the energy qU^ = 18 keV in forward direction 
by a blind (a) on ground potential. 



U -18 kV- 




£16 mm 



Figure 5. The photoelectric electron source (E-GUN). a) ground blind (diameter of 
bore=16 mm), b) gold plated quartz tip at [7^ = —18 kV (yellow), c) metal housing 
(grey), d) insulator, e) sapphire window, f) deuterium UV lamp. The inset shows 
the metal housing which includes the gold plated quartz tip. The ground blind was 
removed here. 



The electron source is mounted on a manipulator which allowed us to move it on a 
sphere up to a ±19° into horizontal and vertical direction (cf. Fig. H]). This corresponds 
to a motion for the gold plated tip on a radius of 1.06 m around the point with the 
pre-spectrometer coordinates z = —2.4 m and r = m, 0.25 m behind the center of the 
source magnet (cf. Fig. H]). The intensity stability of the electron source was measured 
to better than 0.2 % per hour. Fig. [5] shows the cross-section of the E-GUN. The 
reader should keep in mind that the electric field near the gold plated tip (b) is strong 
and there is a ground blind (a) mounted in front of the tip. Together with the electric 
and magnetic fields inside the E-GUN and the detector, these two components play an 
important role in our data analysis (Section EJ. 

The pre-spectrometer detector is a quadratically segmented silicon PIN diode with 
64 pixels of equal size and properties. It has an overall sensitive area of 16 cm 2 [14J. 
It is a predecessor of the final KATRIN detector and was manufactured with the same 
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processing techniques. For simulations of the detector response, the energy resolution 
^•^FWHM an d the dead layer thickness A have to be known. The detector system 
exhibited a measured average A-&pWHM ~ 3.5 keV. The thickness of the dead layer 
was determined as A = 119 ± 3 nm and A = 109 ± 3 nm by using two independent 
experimental techniques [15] . The detector is located at z — 2.3 m (15 cm behind the 
center plane of the detector magnet) at B = 3.4 T and at ground potential (Fig. 0J. 
If the detector is centered on the PS axis, its area corresponds to 28.5 % of KATRIN's 
magnetic flux tube. In our measurements, the detector was adjusted laterally so that 
only a single pixel was hit and the data analysis was made with this single pixel. The 
energy calibration of the pixel was used to select events in the region of interest from 
15 keV to 21 keV. 

For our measurements, the inner electrodes (E) and (F) (cf. Fig. H]) were put on 
a voltage of ^electrode — ^tank = — 500 ±0.1 V with respect to the PS tank using a 
voltage supply (Canberra 3101/2) mounted inside a rack on tank potential. The tank 
voltage itself (U^-^ between -0.5 kV and -17.5 kV) was supplied by another voltage 
source (FUG HCN 140M-35000). The gold plated tip of the E-GUN was supplied with 
a constant voltage of U^ = —18 kV by a high voltage supply (FUG HCN 35-35000). 
The accuracy of the voltage difference determination between the tank and the gold 
plated tip was better than 10 V. 

4. Measurements 

We performed six measurement series with E = 18 keV electrons (Tab. [1]). The PS 
solenoids were set to the KATRIN design value of -B so ] = 4.5 T (Figs. [T]and|l]) and half 
this value -B go ] = 2.3 T. As the cyclotron length grows with 1/B (Eq. [T|), deviations 
from the ideal transmission properties described in Section |2] are more probable for the 
decreased B-field B ^ = 2.3 T. The electrons pass the PS on a radius r(z), which 
encloses a constant magnetic flux 0enc (in a homogeneous B-field one has 0enc ~ 
it ■ r 2 (z) ■ B(z)). For each B-field, the PS transmission was measured for three different 
angular positions (a e {0°, 15°, 19°}) of the E-GUN (cf. Fig. 0J. With the given angles 
a, electrons pass the central plane of the PS at the radii r ce n £ {0 cm, 41.4 cm, 52.3 cm}. 
At the KATRIN design value -B go j = 4.5 T, the corresponding fraction / = ^enc/^KAT 
of the enclosed KATRIN flux tube KAT = 191 T • cm 2 is / G {0 %, 45 %, 72 %}. 

In each measurement series, the PS filter energy qU-p^ was stepped repeatedly from 
1 keV towards 18 keV and back to 1 keV, using identical time intervals. Combining the 
counts from the ramp up and the ramp down eliminates a possible linear drift in the 
emission rate of the E-GUN. Each detector run at a constant potential lasted for about 
R = 48 s. The procedure was repeated up to five times (see 'scans' in Tab. H}, resulting 
in an overall measurement time of e.g. 478 s ~ scans • 2 • R for the measurement with 
-B go j = 4.5 T and r C en = cm. From the runs with the identical voltages, r C en and 
B-field settings, electron events were summed and divided by the overall measurement 
time to obtain an average electron rate. The detector was adjusted laterally before the 
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Table 1. Summary of measurements. 
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start of the measurements so that only a single detector pixel was hit. Only events 
within the region of interest 15 keV < E < 21 keV and from this single detector pixel 
were counted. The comparison of our six measurement series at keV surplus energies 
with our simulations are shown in Figs. |6] and [71 Except for a single measurement 
(-B go l = 4.5 T, r C en = cm), the detector rates at positive surplus energies always 
decrease with growing surplus energy E — qU-pQ- The explanation of this observation is 
given in the section |6j 

5. Simulation Tools 

The analysis of our experimental data is done by comparing it to computer simulations. 
The main components of our simulations are: (i), electric and magnetic field 
computations; (ii), electron tracking in vacuum; (iii), electron scattering with H2 
molecules; (iv), electron tracking in silicon. For the electric field calculations 
(axisymmetric and three-dimensional, with wires) we used the boundary element method 
[j~6| ITT] . In order to speed up the simulation in axisymmetric regions, we employed the 
zonal harmonic expansion method [TTJ. This turned out to be useful also for the E- 
GUN - pre-spectrometer geometry, which is not axially symmetric but consists of two 
separated, locally axisymmetric regions (E-GUN region and pre-spectrometer region). 
The zonal harmonic expansion method was also used for magnetic field computations 
[j~2] . For the electron tracking in vacuum, the exact relativistic equation of motion of the 
electron with Lorentz force was employed [TTJ, using an explicit 8th order Runge-Kutta 
method to solve the ordinary differential equations |18j . The electron-Id^ scattering 
code contains total and differential cross sections p^l f^Dl I21T [22j |23| f2^| l25j f26l f2T] 
and Monte Carlo generation algorithms for elastic, electronic excitation and ionization 
collisions of electrons with H 2 molecules [TTJ. Electron detection, electron energies 
deposited in the sensitive volume of the silicon detector, the detector dead layer, and 
electron backscattering at the detector are modeled by a Monte Carlo C++ code 
(KESS: KATRIN Electron Scattering in Silicon), which is based on detailed studies 
[28| 1291 [301 l3T] and agrees well with experimental data [32] . 

Our original field calculation, 8th order Runge-Kutta tracking and e-H 2 scattering C 
codes pZTJEGldT] have been rewritten into C++ and integrated into the global KATRIN 
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Figure 6. Measurement (blue circles) and simulation (red squares) for -B go i = 4.5 T. 
The simulation was normalized so that the rate at E sm - = E — qllpc = 0.5 keV 
is 1. The experimental rate was normalized so that its average rate is equal to the 
average of the simulated rate. Only simulation values where a measurement exists 
were considered for this average. The simulation points are shifted to the right by 
0.2 keV to make the points distinguishable. The plot shows that our simulation and 
our measurements are compatible at the percent level. 
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Figure 7. Measurement (blue circles) and simulation (red squares) for B go j = 2.3 T. 
The simulation was normalized so that the rate at Esur = E — gfpg = 0.5 keV 
is 1. The experimental rate was normalized so that its average rate is equal to the 
average of the simulated rate. Only simulation values where a measurement exists were 
considered for this average. The simulation points are shifted to the right by 0.2 keV 
to make the points distinguishable. At a = 19° / rcen = 52.3 cm we measured at 
fewer filter potentials. The plot shows that our simulation and our measurements are 
compatible at the percent level. 
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C++ simulation framework 'Kassiopeia' [23]. KESS [31] has also been integrated into 
Kassiopeia. We have used both the original C codes and the new Kassiopeia C++ code 
for the simulations of our paper. 

6. Simulation and Analysis 

Ignoring the influence of the E-GUN and the detector, the PS transmission at keV 
surplus energies in our codes [TH [121 E] is always 100%: In order to show this, 
we started 130 electrons with uniformly distributed polar- and azimuthal angles for 
B sol G {2.3 T,4.5 T}, t/ electrode = ^ tank -0.5 kV G {-1.5 kV, -2.5 kV, .., -17.5 kV} 
and ?" s tart *= (^ mm > 26 mm, 37 mm} in the entry side magnet (z = —2.15 m) of the 
PS. For every simulated electron, the exit condition (z > 2.15 m and r < 37 mm) 
was reached. In the following we show that the measured decrease of detector rate at 
positive surplus energy E — qU-p^ is not caused by a loss of transmission in the PS, but 
by the special setup (Fig. H]) used for the measurements. The transmission probability 
of 100% can only be confirmed indirectly with an uncertainty at the percent level via 
the comparison with a simulation. The essential ingredients to explain our data are: 

(i) A constant PS transmission probability of 100% 

(ii) Electron backscattering at the detector 
(iii) B x V-B-drift of electrons proportional to -E^in ( see ©) 
(iv) Reflection of electrons in the electric field of the E-GUN 

(v) The loss of backscattered electrons hitting the ground blind in front of the E-GUN 
(Figs. [5] and ED 

Electrons impinging on the detector with an incident energy E^ nc = 18 keV and an 
incident angle 0- inc = 0° have a probability of about 20% to be backscattered, and higher 
incident angles further increase this probability [351 ES] . Most of the backscattered elec- 
trons have lost energy in the detector and are again reflected by the filter potential 
gf/pg in the PS or by the magnetic mirror effect towards the detector. Backscattered 
electrons having deposited less than E — qU-p^ in the detector retain enough energy to 
pass the filter potential qUpg of the PS in backward (towards E-GUN) direction. They 
can enter the E-GUN through the opening in the ground blind (cf. Fig. |5]) and can be 
electrically reflected towards the detector again. This process continues until all energy 
is finally deposited inside the detector. Even at high surplus energy E — qUp^, the 
count rate at the detector is constant. Flight times in the PS are of the order of 10 ns 
which is far below the /xs shaping time of the DAQ. Therefore, only electrons with large 
energy losses in the deadlayer will deposit an energy lower than the region of interest in 
the sensitive volume. Since all electrons hitting the detector have the same energy and 
same angular distribution for all settings of ^ e i ec trode = ^tank - ^-5 kV, the count rate 
does not depend on this voltage setting. This explains the measurement at -B go j = 4.5 T 
and r cen = cm (a = 0°), where no loss in count rate is observed. 
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The description above is also valid for the off-axis E-GUN settings with r cen G 
{41.4 cm, 52.3 cm} (a G {15°, 19°}). 

The axial rotation (J7J can cause the electron to eventually hit the ground blind 
(Fig. IS]), depending on B, the electron energies E± and E\\ and the total path length. 
The axial rotation is only dependent on the energy and has the same sense for a for- 
ward (towards detector) and backward (towards E-GUN) pass of the PS and is therefore 
adding up at each pass of the PS. The path for the electrons in the pre-spectrometer is 
elongated by reflections at the detector, at the field of the E-GUN, at magnetic fields 
and at the spectrometer potential and can thus be multiples of the spectrometer length. 
The higher the surplus energy, the higher the probability for a backscattered electron to 
overcome the spectrometer potential after energy deposits in the detector. The higher 
the surplus energy, the larger the axial rotation which guides electrons onto the ground 
blind. Therefore, the count rate in the region of interest decreases with higher surplus 
energy -Esur = Eq — ^PS' With this, all measurements at -B go ] = 4.5 T can be ex- 
plained. 

For -B go } = 2.3 T a loss in count rate is observed for all E-GUN settings, including 
the one at r cen = cm (a = 0°). It is not possible to explain this effect with 
backscattering and the B x V|S| drift alone. The electric field gradient in the E- 
GUN is large compared to the PS, since the potential difference of 18 kV is applied 
across a distance of only a few cm. Together with the 50% lower magnetic field in the 
center of the magnets, this can lead to a non-adiabatic transport in the E-GUN region. 
At S gol = 2.3 T, the E-GUN is located at a magnetic field of B « 0.02 T. Thus, a 
backscattered electron entering the E-GUN through the ground blind has a probability 
to change its angle towards the magnetic field line non-adiabatically. Depending on 
the new angle and the electron energy, it can be trapped between the E-GUN and the 
closest magnet or the spectrometer potential. Thus, a loss in count rate in the region 
of interest will also be observed for measurements with r cen (a = 0°). This assumption 
is therefore able to explain the measurement for _B go j = 2.3 T and r ce n (a = 0°). 

In our simulations, the electrons were started with a uniform random kinetic 
energy of < E < 2eV, and were uniformly distributed on a disc with diameter 
d = 1mm in front of the actual gold tip. An angular distribution, 9 = arcsin(i?i) 
from [37] and ip = 2ttR 2 with uniformly distributed random numbers Ri, R 2 G [0, 1[ 
was used. The time an electron travels between two subsequent detector hits is more 
than two magnitudes smaller than the DAQ shaping time. This means, subsequent 
hits are analyzed by the DAQ as one hit. Therefore, each energy deposition in the 
sensitive detector volume per electron was summed up, even for electrons with multiple 
detector entries. As in the experimental data analysis, all electrons with energies 
15 < E < 21 keV were counted. 
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Figure 8. The azimuthal drift in the pre-spectrometer. Left: Geometry of the E- 
GUN. Right: The azimuthal drift seen from the detector side. The electron starts in 
the centre of the bore in the ground blind. Towards the axis of the pre-spectrometer, 
the B-field is stronger and the curvature of the electron trajectory is larger, resulting 
in an overall B x VB drift around the pre-spectrometer axis ([7]). If the drift is large 
enough, the backscattered electron hits the ground blind on its way back from the 
detector. 



Possible exit conditions of the simulation were: 

• Electron hits E-GUN blind (cf. Fig. |5]and Fig. E>. 

• Electron energy lower than 100 eV. It is no longer able to pass the dead layer of the 
detector. 

• Electron was reflected more than 20 times in the PS (trapping). 

Figures |6] and [7] show the normalized simulation and experimental results. They 
show that our simulation can explain the experimental rates at the percent level. 
The agreement between experiment and simulation implies that the PS has 100% 
transmission probability for electrons with keV surplus energy as long as B SQ \ is not 
smaller than 2.3 T (half the KATRIN design value pQ). 



7. Operating the PS at reduced zero filter energy 

The previous section proves experimentally that the PS can be operated at reduced or 
even zero filter energy without essential transmission losses, as long as B ' j is not smaller 
than 2.3 T. In this section we discuss more generally, with the help of simulations, the 
operation of PS with zero potential. We show first that, with f/pg = , the background 
due to positive ions created by beta electrons is expected to be far below the KATRIN 
requirement of b = 10 mHz, and second that the electron motion through the PS and 
MS is fully adiabatic. 

Lowering the filter energy qU-pg of the PS, the depth of the natural Penning trap, 
created by the PS and MS retarding potentials f/pg and Uy^a and the magnetic field 
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£> so l, between the PS and MS decreases, and the corresponding background level is 
expected to decrease, too. On the other hand, there is another background component 
that increases with decreasing qU-pg: more /3-electrons reach the entrance of the main 
spectrometer. These electrons are not able to produce direct background, but they create 
positive ions through ionizing collisions with the residual gas molecules. These ions can 
fly deep into the main spectrometer, and they can produce low energy electrons there, 
either in the residual gas or at the inner surface of the main spectrometer electrodes. 
Some of these electrons can hit the detector, and so we obtain background events. We 
will present below a quantitative estimate for the maximal value of this background 
component, in case of zero PS potential. 

Inside the WGTS, 10 11 electrons are produced through tritium beta decays each 
second. About 20 % of them move through the transport system and reach the 
entrance of the main spectrometer in the case of qU-pg = keV. Along their way, 
many of these electrons have ionizing collisions with residual gas molecules and thus 
produce positive ions. In the absence of any hindrances, all these ions enter the main 
spectrometer. In order to reduce the background rate of these ions, we plan to use an 
ion-blocking electrode near the center of the main spectrometer entrance magnet. Since 
the positive ions created by ionizing collisions have small kinetic energy (below few 
eV), an accordingly small potential barrier created by an ion-blocking electrode would 
already prevent these ions from entering the MS and thus from producing background. 

Nevertheless, the ions created in the region between the ion-blocking electrode and 
the high- voltage area of the MS are able to enter the MS. We have computed the ion 
creation rate in this region by detailed trajectory simulations. First, we generated 
electrons at the center of the MS entrance magnet by using the tritium /3-decay 
Fermi spectrum and an isotropic angular distribution. These electrons were tracked 
until their reflection at the MS filter potential and then back to their starting point. 
The ionization probability of an electron was computed by summing the differential 
ionization probabilities dV loa = <Ji on (Ekm) -n-dl, where cr iori (-£7kin) denotes the ionization 
cross section of electrons with the residual gas molecules, as function of the electron 
kinetic energy E^ in , n is the number density of the residual gas, and dl the differential 
pathlength. In our simulations, we assumed molecular hydrogen for the residual gas, 
and we used the e — H 2 ionization cross section formulas of [20j E2] (they are in good 
agreement with measured e — H 2 cross section values). Assuming p = 10~ n mbar 
pressure and room temperature, the number density is n = 2.4 ■ 10 11 m" 3 . According 
to our calculations (simulation of 1000 electron tracks), the average ionization cross 
section is a lon = 10~ 21 m 2 , the average electron pathlength is I = 1.5 m, and the average 
ionization probability of an electron is Vi on = 3.6 • 10 -10 . Using the 2 • 10 10 s _1 /3-electron 
intensity, we obtain a positive ion creation rate of roughly iV + rj 7 s _1 in the region 
between the ion-blocking electrode and the high potential domain of the MS. 

If an electron is scattered towards large polar angles 9 (remember E± = ^j n -sin 2 9), 
it can become trapped in a hybrid trap near the entrance of the MS: if this electron 
moves towards the entrance magnet of the MS coming from inside the MS, 9 increases 
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adiabatically until 9 = 90° is reached, and the electron starts to move towards the MS 
again. The MS magnet thus establishes a magnetic mirror for these electrons. Inside 
the MS, the electron is electrically reflected by the MS filter potential generated by 
the MS electrode system [38]. In order to compute the ion creation rate due to these 
trapped electrons, we simulated 10 8 electrons. We used our custom C codes [TO HU [12] 
to compute the electromagnetic fields, the trajectories and the e — H 2 scattering. The 
result of these simulations is the following: the average trapping probability of the beta 
electrons in the hybrid trap is 3 • 10 -11 , and the number of ions created by a trapped 
electron is smaller than 5 (the trapped electrons can leave the trap by scattering and 
by energy loss due to synchrotron radiation). From these numbers and from the above 
beta intensity we get an ionization rate that is smaller than 3 s" 1 . Therefore, the ion 
creation rate due to these trapped electrons is smaller than due to the free (non-trapped) 
electrons. 

The ions with roughly N + ps 10 s" 1 rate enter the MS, they will be accelerated to 
about 18.57 keV kinetic energy, and due to the small (few Gauss) magnetic field, their 
motion inside the MS is completely non-adiabatic: they move on a straight line, until 
they hit the spectrometer tank. During this motion, they can suffer ionizing collisions 
with the residual gas. The ionizing collision cross section of H + and H^ ions of 18 
keV kinetic energy with H 2 molecules is about a^ n = 2 • 10~ 20 m 2 [39]. Assuming 
l + = 20 m pathlength for the positive ions inside the main spectrometer tank, the 
secondary electron creation rate due to ionizing collisions of the positive ions with H 2 
molecules is N e = N + a^ Q l + n = 10~ 6 s _1 , corresponding to 0.001 mHz background 
level. This background increases quadratically with the residual gas pressure, so with 
p = 10~ 10 mbar the background rate would be 0.1 mHz. 

Another background possibility is the following: the ions hit the MS tank with 
high velocity, and these impact events are connected with secondary electron emission 
from the surface. One ion can eject more than 1 electron; let us assume that this 
multiplication number is 10. Then we get a secondary electron emission rate of 100 s _1 
from the tank surface. This is several orders of magnitude smaller than we expect from 
cosmic ray muons and environmental radioactivity. Thanks to the magnetic shielding 
of the approximately axisymmetric magnetic field and to the electric shielding of the 
wire electrode, only a very small proportion of these electrons is expected to reach the 
detector; extrapolating the experimental data of the Mainz neutrino mass spectrometer, 
this proportion could be about 10~ 7 . With this suppression factor, we get 0.01 mHz 
background level from these electrons. 

For both scenarios we obtain a background level caused by positive ions which is 
several orders of magnitude smaller than the b = 10 mHz background value that would 
be acceptable for the KATRIN experiment [lj. Therefore, our simulations show that 
the PS could be used with zero potential, without any significant background increase 
due to the positive ions produced by the beta electrons. 

With zero or small PS filter energy, the signal electrons (Eq = 18.57 keV) have 
a high surplus energy E SUT = Eq — gf/pg inside the PS. Due to the relatively small 
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magnetic field in the middle of the PS (-B ana pg = 0.016 T), it could happen, in 
principle, that the motion of these electrons is not adiabatic. We have checked the 
adiabaticity behaviour of the electrons with detailed trajectory simulations. For this 
purpose, we started the electrons in the KATRIN source (WGTS) at various points 
with various polar and azimuthal direction angles, and we tracked them as far as 
the main spectrometer analyzing plane. We defined the starting kinetic energy of the 
electrons with the following procedure: first, using the starting point, direction vector 
and a first estimate for the transmission energy, we computed the guiding center point 
corresponding to the starting point. Then, we simulated the magnetic field line from 
the guiding center point until the MS analyzing plane. Using the electric potential 
and magnetic field values at the two endpoints of this field line, and the starting polar 
angle, we calculated the adiabatic transmission energy: in the adiabatic approximation, 
the electron has zero longitudinal energy E\\ in the analyzing plane (cf. Fig. [2]) if it 
starts with this energy (electrons starting below or above this energy are reflected or 
transmitted, respectively). We defined the starting kinetic energy of our electrons as the 
above adiabatic transmission energy plus a small surplus energy (E SUT = E — g£/pg = 3 
meV). If the electron motion is fully adiabatic, E\\ has to be precisely equal to E SUT in 
the analyzing plane. The main result of our simulations is the following: both for a PS 
filter energy of qUp^ = keV and qUpo = 18.3 keV, and for all starting parameters, 
using the standard KATRIN magnetic field values (3.6 T in WGTS, etc.), \E\\ — E sur \ 
(computed by exact tracking) is in the MS analyzing plane on the average 0.2 meV, 
which is four orders of magnitude smaller than the resolution A-E^g = 0.93 eV of the 
KATRIN MS. Therefore, we can say that the motion of signal electrons in the KATRIN 
system is practically adiabatic, even with Upc = kV; deviations from adiabaticity have 
a negligible effect to the KATRIN transmission function. Note that that the electron 
motion is approximately adiabatic even if the magnetic field in the whole KATRIN 
system is half of its standard design value; in this case, \E» — E snr \ is in the MS analyzing 
plane on the average 0.6 meV. 

The local behaviour of the quantity 7/x was considered as adiabatic invariant in 
Section [2j There, the difference between the electron energy E and the filter energy 
qUf was assumed to be small (some tens of eV). As one can see in Fig. [U 7// has 
an oscillational behavior, due to the superposition of the cyclotron motion and the 
azimuthal magnetron motion; the oscillation period in Fig. [9] is equal to the electron 
cyclotron period. The amplitude of the oscillation depends on the electron surplus 
energy E snT = E — qUpv: With E snT w 0.3 keV in the PS, the relative fluctuation 
Sjfjb/jf/, inside the PS is order of 10 -2 (for smaller starting polar angle 9 the fluctuation 
is somewhat larger), but with E sur of a few keV this fluctuation is much larger, order 
of 1 (Fig. ED. 

How is it then possible that 7/x regains its starting value, in spite of its large 
oscillations in the small magnetic field region? The explanation is the following 
[10~| HU 02]: the real adiabatic invariant J a j, which is constant in the adiabatic 
approximation throughout the whole trajectory, is not 7/z, but a complicated function 
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Figure 9. Oscillations of yn and 
regain of original value: B so \ = 
4.5 T, E/ ps = -1.5 kV, E = 
18 keV, z start = -2.15 m, <^ start = 
60° (azimuthal angle), s tart = 
72° (polar angle) and r start = 
26 mm. The PS solenoids are at 
z = -2.15 m and z = +2.15 m. The 
radius of the KATRIN flux tube in 
solenoids is 37 mm at 5 g0 j = 4.5 T. 
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Figure 10. Evolution of 7/1 for a 
trapped electron: B g0 ] = 0.57 T, 
t/pg = -1.5 kV, E = 18 keV, 

2 start = -2- 15 m > ^start = 30 ° 
(azimuthal angle), # s tart = ^2° 
(polar angle) and fgtart = 26 mm. 
The PS solenoids are at z = 
-2.15 m and z = +2.15 m. The 
radius of the KATRIN flux tube in 
solenoids is 37 mm at B SQ \ = 4.5 T. 




Figure 11. Evolution of the polar 
angle for a transmitted electron. 
Parameters as in Fig. [9] 



Figure 12. Evolution of the polar 
angle 9 for a trapped electron. 
Parameters as in Fig. QUI 



of the higher field derivatives. Inside homogeneous field regions the field derivatives 
are small, and there Igj ~ 7/i. When an electron moves from one homogeneous field 
region to another, and if the adiabatic invariant J„j is constant, 7/1 can regain the 
starting value with high accuracy, although between these two regions, where the field 
gradients are large, 7/1 can have large oscillations. Using a human analogy, we can say 
that the electron seems to have 'memory', remembering its initial value of 7/i jl2]. If 
the electron motion is not adiabatic (in case of high energy or small magnetic field), 7 ac j 
is not constant, so in this case the electron will not gain back its original 7/1 value in 
the second homogeneous field region; then, the electron has no memory. 

In order to illustrate this phenomenon, we simulated electrons starting in the entry- 
side magnet of the PS. At the KATRIN field 5 gol = 4.5 T (Fig. [l]and|4]), electrons 
are transmitted through the PS, and 7/i regains its original value (Fig. [9]). If we choose 
a deliberately low field B ' j = 0.57 T in simulations, electrons can be magnetically 
trapped and 7/i becomes chaotic (Fig. [10]). If 7yU regains its original value, the angle 6 
is approximately determined by Eq. [3] Let us consider an electron which enters the PS 
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in the entry-side magnet with 9 < 90° and for which the polar angle 9 is determined by 
Eq. [3j In this case, the electron will never acquire 9 = 90° anywhere in the PS. Thus, 
it cannot be magnetically reflected (Fig. [IT]) . On the other hand, if Eq. [3] is violated as 
shown in Fig. [10J the reflection angle 9 = 90° can be reached over and over again (Fig. 
IT21) . the electron is magnetically trapped inside the PS. It is not transmitted, and so it 
is not detected. 

We mention that, at the edge of outer field lines of the flux tube, with £? so i=4.5 
T and zero PS potential, the 18 keV electrons make inside the PS about 2 degrees 
magnetron motional rotation around the beam axis. This has to be taken into 
account for precise imaging investigations of the KATRIN experiment (see Ref. jl] 
for experimental examples of much larger magnetron motional rotation). 

The background as a function of the PS potential will be investigated experimentally 
when the whole KATRIN system is finished. Similarly, the MS transmission function and 
thus the adiabaticity of the electrons in PS and MS can be experimentally investigated 
rather precisely, as function of PS potential, if PS and MS are connected together (by 
shooting E-GUN electrons through them). The experiments presented in our paper, 
using only the PS, are sensitive only to large deviations from adiabaticity (small non- 
adiabaticity effects do not cause any transmission losses in the case of large surplus 
energies) . 

8. Conclusion and Outlook 

Our investigations show that the PS filter energy <?£/pg can be reduced by several 
keV without any loss of transmission, making it possible to diminish the Penning 
trap between PS and MS. The actual value of the PS filter energy qllpg to minimize 
KATRIN's background has to be determined experimentally with the full KATRIN 
setup. It is possible that the final KATRIN setup will operate with a mixture of 
reduced PS filter energy and other, active measures removing stored electrons from 
the Penning trap between PS and MS. Sweeping a wire through the trapping volume 
[5] has proven to be an efficient means to empty the trap. It will take about £ w j re = 1 s 
for the wire to sweep across the magnetic flux tube imaged on the detector. In order to 
scan the tritium /3-spectrum, the MS retarding voltage Uy^a will be changes every few 
minutes. The sweeps will be performed during these voltage changes to avoid a loss of 
measurement time. Yet, the corresponding 10 6 motion cycles are a very large number 
for an UHV compatible device. 

For KATRIN, the transmission function of the PS has to be known with permille 
accuracy pQ. Our investigations indicate that this precision is only possible if the 
influence of multiple electron reflection at the electron source is suppressed. This can 
be achieved with a) a stable, pulsed electron source at the entry of the PS, b) the PS 
and MS in their final tandem configuration (cf. Fig. [1]) and c) KATRIN's final detector 
having a time resolution of about 100 ns. 

A novel, angular selective pulsed UV laser photoelectron source, which can produce 



The KATRIN Pre- Spectrometer at reduced Filter Energy 21 

pulses as short as 40 ns with a repetition rate of up to 10 kHz, is currently being built 
for KATRIN on the basis of [43j HU H5]. In the experiment we propose, the MS will be 
operated so that the electrons retain only a few eV surplus energy in the MS and are 
therefore guaranteed to be transmitted if the MS performs as it should. The filter energy 
of the PS will be varied, just as in our experiment. In this configuration, multiple (up 
to about 14) reflections between the detector and the main-spectrometer potential will 
cease after less than 20 /xs. The fraction of backscattered electrons, which retain enough 
energy to pass the MS filter potential in backwards (towards electron source) direction 
and could possibly get lost at the electron source, is negligible in this configuration. 
Therefore, the analysis of this experiment will be much simpler than in our case. 
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